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', Abstract: We discuss the impact of the planned upgrades of the HERA collider 

on the study of open heavy flavour and quarkonium production. New experimental 
Oh" techniques in charm physics are presented. 
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1 Introduction 

The ep collider HERA offers new opportunities to study the production mechanism of heavy 
quarks and to test the predictions of the underlying theory. In the following, we will deal with 
^ \ the processes 

O ■ _ 

\0 • e + p — > Q + Q + X, (1) 

i + P — Q + Q + x, (2) 

7 + p — (qq) + x, (3) 

where \QQj indicates a generic heavy quarkonium state, like J/if). The dominant contribution 

to the QQ cross section is due to those events in which the virtuality of the photon exchanged 
between the electron and the proton is very small. In this case, the electron can be considered to 
be equivalent to a beam of on-shell photons, whose distribution in energy (Weizsacker- Williams 
function [![]) can be calculated in QED. The underlying production mechanism is therefore a 
photoproduction one (eqs. (@) and (|3|)), which has been studied extensively in fixed-target ex- 
periments. At HERA, the available center-of-mass energy is about one order of magnitude 
larger than at fixed-target facilities; this energy regime is totally unexplored in photoproduc- 
tion, and several new features have to be taken into proper account. In particular, the large 
contribution of the resolved photon component introduces in the theoretical predictions a source 
of uncertainty which is totally negligible at fixed-target energies. 

A complementary way of studying heavy flavour production at HERA is to retain only those 
events characterized by a large photon virtuality (eq. ([!])). Although the total rates are much 
smaller than the photoproduction ones, the resolved component is completely eliminated and 
more reliable theoretical and experimental results can be obtained. Also, the dependence of 
the data upon the photon virtuality can be used as a further test of QCD predictions. 

1 To appear in the proceedings of the workshop Future Physics at HERA, eds. G. Ingelman, A. De Roeck 
and R. Klanner, DESY, Hamburg, 1996. 



2 Theoretical summary 



Fully exclusive, next-to-leading order perturbative QCD calculations @, 0, @] are now available 
for the processes of eqs. (|l|) and (|2]). In the framework of the factorization theorem of ref. ||, 
a next-to-leading order calculation for the direct color-singlet photoproduction of quarkonium 
has been presented in ref. 0. Other contributions to the production of quarkonium, like the 
resolved photoproduction and the direct color-octet photoproduction, which are potentially 
relevant in kinematical regions accessible at HERA, have been considered at the leading order 
in ref. @. 

In spite of these substantial progresses in understanding the production mechanism, the 
(fixed-order perturbative) results may become unreliable in certain kinematical regions, due 
to the appearance of potentially large logarithms which spoil the convergence of the perturba- 
tive expansion. In this case, a resummation to all orders of these large logarithms has to be 
performed. 

When the available center-of-mass energy \fS gets large, the effective expansion parameter 
of the perturbative series becomes a s \og(S/m 2 Q). The problem of resumming these terms 
(small-x effects) has been tackled by several authors ||, mainly in the context of b production 
in hadronic collisions. Specific studies for HERA |J lead to the conclusion that the total 
photoproduction rates can be increased by the resummation of the \og(S/m^)) terms by a 
factor smaller than or equal to 40% (45%) with respect to the next-to-leading order prediction 
for the direct (resolved) contribution in charm production. In the case of bottom, the effect is 
much milder, being always smaller than 10%. 

The transverse momentum distribution is in principle affected by the presence of \og(p T / mo) 
terms. These logarithms can be resummed by observing that, at high p T , the heavy-quark mass 
is negligible, and by using perturbative fragmentation functions [ 10| . Remarkably enough, the 



fixed-order and the resummed results agree in a very wide range in p T for charm production 
(in the HERA energy range, the contribution of the resummation of this kind of logarithms 
is expected to be negligible for bottom production); the effect of the resummation might be 
visible only for p T larger than 50 GeV. 

Finally, multiple soft gluon emission makes the perturbative expansion unreliable close to 
the threshold or to the borders of the phase space, like for example the regions p^9 ~ and 
/\(j)QQ ~ vr. This problem has not been directly dealt with during the workshop, but a lot of 
theoretical work has recently been performed ||11||. The resummation of soft gluons would be of 



great help in order to have a reliable prediction for J/ijj production at p T ~ 0; also, threshold 
effects are important for b production at HERA-B. 



2.1 Quarkonium Photoproduction 

Quarkonium photoproduction has been dealt with in ref. [T2J. The benchmark process is in 
this case the production of J/ijj. With suitable cuts, typically p T > 1 GeV and z < 0.9, where 

k(p) ■ k{ j/ji) 

Hp) ■ Hi) ' 1 ) 

the dominant contribution to the cross section is due to the direct color-singlet production of 
J/i[>] this cross section has been calculated at next-to-leading order in QCD ||. The radiative 



corrections turn out to be sizeable in the HERA energy range, being of the same order of the 
leading-order contribution. The total rate is of the order of tens of nb, and a first comparison 
with data |13] has already been performed. Although there is a good agreement as far as the 
shape of (z and p T ) distributions is concerned, the theoretical predictions appear to slightly 
undershoot the data for the total cross section; this can be adjusted by properly tuning the 
input parameters of the calculation, since the normalization is affected by large uncertainties. 
With a larger integrated luminosity, a more significant comparison will be possible, eventually 
helping in constraining the input parameters. 

The luminosity upgrade is mandatory in order to assess the importance of the J j ip produc- 
tion mechanisms neglected so far, namely the direct color-octet contribution and the resolved 
contribution. The former has been considered in refs. [ff], |12| at the leading order in QCD. It 
has to be stressed that this contribution can not be calculated from first principles, and relies 
upon fits to Tevatron data on color-octet J/ip production. The comparison with HERA data 
can therefore be regarded as a test of the underlying factorization picture [El]. Using the result 
of ref. J7|, the color-octet contribution appears to be dominant over the color-singlet one in 
the large- z region; the first HERA data do not support this behaviour. Nevertheless, since 
both the data and the theoretical calculations need some refinement, no definite conclusion can 
be drawn at present. The resolved J ftp production is expected to be dominant at very low- z 
values, being enhanced by the color-octet terms. Since the total rate is (9(1 nb) in that region, 
large integrated luminosity is needed to study this problem. 

To further investigate the interplay between color-singlet and color-octet contributions, 



other possibilities were taken into account in ref. ||12|| . For example, the photoproduction cross 
section of x should be dominated by color-octet terms; the total rate is however suppressed 
by two orders of magnitude with respect to the one of J/ip, and therefore large luminosity is 
required. The + 7 production has a very distinctive signature to color-octet contribution in 
the large- z region 0, but the cross section is very small ((9(10 pb)). Finally, large luminosity 
is necessary also in order to investigate J/ip production for p T > 10 GeV, where the dominant 
contribution is expected to be due to charm fragmentation (the total rate is at most (9(1 pb)). 



2.2 Open Heavy Flavour Production 

Several aspects of the deep-inelastic production of heavy flavours at HERA have been studied 
during the workshop. The work of ref. basically falls into three parts. To begin with, 
the authors present a study of the heavy quark inclusive structure function F 2 (x , Q 2 , m 2 ) , its 
sensitivity to the input parameters of the calculation, and the size of the QCD corrections. The 
full next-to-leading order theoretical prediction |3|, |j] suffers from a fairly little uncertainty, 
and the radiative corrections are not too large. The structure function is on the other hand 
significantly sensitive to the small-x behaviour of the gluon density in the proton. Therefore, 
it is concluded that F 2 in charm production is an excellent probe to infer the gluon density in 
the proton at small x. 

Next, the exclusive properties of the final state are investigated, since they allow a deeper 
understanding of the dynamics of the production process with respect to the inclusive structure 
function. The next-to-leading order corrections to Fk, k = 2, L, in a fully differential form have 
been recently calculated [|J using the subtraction method. The results were incorporated in 
a Monte-Carlo style program which allows one to study correlations in the laboratory frame. 
The authors of ref. |15j present some single-inclusive distributions and compare them with 



preliminary ZEUS data [16|. In general, the data follow the shape of the theoretical curves, but 



lie above them. This behaviour is consistent with the recent analysis by HI |L7j that, using a 
single-inclusive quantity, showed that the charm production mechanism is boson-gluon fusion, 
and not a consequence of charm quarks in the sea. Correlations between charm and anticharm 
are also presented. Clearly, all these issues will be further clarified as soon as data with better 
statistics will be available; a large integrated luminosity is highly desirable. 

Finally, the authors of ref. |L5] deal with the fact that, at truly large Q 2 , a charm quark 
should be described as a light quark, i.e. as a constituent parton of the proton, whereas at Q 2 



of order of m 2 it can only be produced through boson-gluon fusion (see ref. [T8| for a detailed 



discussion). It is shown that, surprisingly, already at Q 2 of order 20-30 GeV 2 , the asymptotic 
next-to-leading order formula (large Q 2 ) for the inclusive structure function differs from the 
exact result for a factor of 5% or less, indicating that at these not so large Q 2 values the charm 
quark behaves already very much like a parton. Since, as mentioned before, the boson-gluon 
fusion almost correctly accounts for experimental results on single-inclusive distributions, this 
points out that the scale at which the charm quark can be regarded as a parton actually depends 
upon the physical observables considered. 

Open heavy flavour photoproduction has been considered in ref. With an integrated 
luminosity of 1000 pb~ x , next-to-leading order QCD predicts about 10 9 (10 6 ) charm (bottom) 
particles produced in ep collisions (Weizsacker- Williams approximation). Therefore, taking into 
account the experimental efficiencies, the number of reconstructed particles will be comparable 
to or larger than the one obtainable by the fixed-target hadropro duct ion experiments of the 
new generation at Fermilab. This will provide the HERA experiments with the possibility 
of performing studies of charm and bottom physics at an excellent level of accuracy. In this 
respect, it is extremely interesting to compare the results at HERA with the results of fixed- 
target photoproduction experiments, which have a center-of-mass energy of about one order of 
magnitude smaller. In charm physics, high luminosity will allow to consider exclusive quantities, 
like correlations, which constitute the most stringent test for the underlying theoretical picture. 
One may also adopt a different point of view, namely to look at charm production as a useful 
tool to constrain the input parameters entering the calculations, like for example the quark 
mass. To this end, it is mandatory to have data with large statistics. The capability of 
HERA of producing bottom quarks is also very promising. In particular, the comparison of 



the QCD prediction for the p T spectrum of bottom (which has been shown in ref. [20| to be 
only marginally affected by the uncertainties on the input parameters of the calculation) with 
the data could be of great help in understanding the origin of the discrepancy observed at the 
Tevatron for the same quantity. 



2.3 Production at HERA-B 

The study of the production mechanism of heavy flavours at HERA-B might prove to be 
extremely useful. The data on bottom production at fixed target have low statistics, and have 
been mainly obtained in collisions (only very recently, the first measurement of the total 
bb rate in pN collisions has been presented). A study of distributions would be interesting in 
order to test the QCD description of bottom quark production at low center-of-mass energies. 
In charm physics, the results of fixed-target experiments at CERN and Fermilab still leave 



plenty of open questions (see ref. |21] and references therein) . In particular, new measurements 



may help in understanding the importance and the nature of non-perturbative contributions 



to charm cross section. The hypothesis of the intrinsic charm in the proton could be tested as 
well. 



2.4 Determination of the Gluon Density in the Proton using Charm 
Data 

The possibility has been considered of using charm data to constrain (or to measure) the 
gluon density in the proton, both in the unpolarized and the polarized case. Photoproduction 
of charm is in principle an ideal tool in order to perform this measurement, since the gluon 
density enters the cross section already at the leading order, and in a simpler way with respect to 



hadroproduction processes. It has been shown in ref. [^2] that the use of correlations between 
the charm and the anticharm can give a direct measure of the unpolarized gluon density in 
the proton in a fairly large range in x (10 -3 -v- 10 _1 ), since the contribution of the resolved 
component can be suitably suppressed. Being necessary to reconstruct both the charm and the 
anticharm for this kind of measurements, a very large integrated luminosity is required. The 



measurement of the polarized gluon density in the proton has been considered in refs. E4 



see also ref. [2l|), in view of the possibility for HERA to operate in the polarized mode. The 



conclusion has been reached p3j that charm photoproduction data can be used to constrain 



this quantity, if an integrated luminosity of at least 100 pb 1 will be achieved. 



2.5 Intrinsic Charm 

The possibility of detecting signals due to intrinsic charm in the proton has been discussed 
during the workshop by G. Ingelman. An intrinsic charm contribution has been suggested 
long time ago J26| to explain an excess of data from fixed-target experiments with respect 



to theoretical predictions in the large-x F region. Recently, this problem has been tackled in 
ref. [0, dealing directly with ep collisions at HERA. It turns out that the intrinsic charm 
contribution is non-negligible only in the very forward region, and can not be detected with the 
present experimental configuration; an upgrade in the very forward region would be necessary. 



3 Heavy Quark Production Cross Sections and Recon- 
struction Efficiencies 

In this section we give measured cross sections and expected reconstruction efficiencies for 
charmed mesons. As an example the capabilities of the HI detector [2S| are given. The effect 
of a double layer silicon vertex detector p9f is investigated. 



The charm production cross section has been measured by both collaborations HI and 
ZEUS [BUI pij and is close to one microbarn. The value for photoproduction (Q 2 < 0.01 GeV 2 ) 



amounts to |30[| 



a[ep -> ccX) = 941 ± 160±J|g nb. 

The first error is the statistical and the second the systematic error dominated by the extrapo- 
lation uncertainty to the unmeasured phase space. This rather large cross section makes HERA 
an ideal place for charm physics. 



The b cross section has not been measured yet. An estimate o[ep — > bbX) = 6 ± 2 nb is of 
similar magnitude as on the Z°- or the T-resonance and much lower than at hadron colliders. 

The reconstruction of charm uses the central tracking chambers with an angular coverage 
of —1.5 < 7] < 1.5. The decay chain 

ep -> D*X -> D°7T S -> (K-it + )it s 

yields a signal to background ratio of 1:1 for the mass difference AM = M(K~7t + tt + ) — 
M(K~tt + ) with a cut on the transverse momentum of p T (D*) > 2.5 GeV. A total efficiency 
in photoproduction of trigger, acceptance, branching fraction, and reconstruction of charmed 
events via -D*-tagging of 10 -4 has been achieved (table [I]). This is only a tiny fraction of the 
total cross section. The yield can be enlarged by summing other decay channels of the D meson. 
A similar signal to background ratio can be achieved with the reconstruction of 

D° -> K~tc + (3.83%), K\ + ti- (5.4%), R-n+n+n- (7.5%), K~fj+u (3.23%). 

This is an improvement of a factor of 5 over the single K~ir + channel. 

Besides the branching fractions into a specific decay channel the biggest event loss (order of 
magnitude) comes from the transverse momentum cut of the D* which is necessary to reduce 
the combinatorial background. This can be partially compensated with the help of a vertex 
detector. 



3.1 Charm Reconstruction with the HI Vertex Detector 

The main effect of a vertex detector is the assignment of charged particle tracks to a sec- 
ondary vertex which results in a much improved signal to background ratio S/N, but with a 
corresponding loss of acceptance. 

This can be seen in figure |l]a (from ref. ||33|Q, where this ratio S/N is plotted versus the 
separation d/a of primary and secondary vertex, with d the distance and a its error. No lower 
momentum cut for the D* was applied. The track finding efficiency of the slow pion in the 
D* — > D°7r s decay is non zero only above p T (7r s ) = 120 MeV which implies an indirect transverse 
momentum cutoff for the D* meson of about 1 GeV. 

The signal to background grows rapidly with the vertex separation and reaches values above 
one for d/a > 2 — 2.5. At the same time we observe a loss of the signal by a factor two compared 
to no vertex cut (see figure [L]b). Therefore the net gain is not very large. 

But with the vertex cut even more decay channels are reconstructable. The D° mesons 
can be found directly without the detour via the D*. This gains a factor of more than two. 
Also the D+ decays into K-it + it+ (9.1%), K°ti + it + it- (7%), and K°n + (2.8%) will be possible. 
Therefore the net gain with a vertex detector to fully reconstruct charmed final states is a factor 
5-10. 

A vertex detector is also useful for tagging charm by looking for secondary vertices inside 
jets. A study |33] has been performed where two or more jets with E ] T et > 6 GeV, \r]^ et \ < 2 , 
Q 2 <4 GeV 2 , and 135 GeV < W yp < 270 GeV were selected. The efficiency using very loose 
cuts (at least one track with an impact parameter displaced by d/a >2) is 37% for charm and 
17% for bottom with S/N=0.9 and also here the gain from a vertex detector is significant (order 
of magnitude). 

Tables |T] and |^ summarize the efficiencies. For an integrated luminosity of 300 pb~ x we 
expect therefore of order 10 6 fully reconstructed D° and mesons with a S/N>1. 




Figure 1: a): Signal to background ratio S/N of the D* signal versus the vertex separation d 
normalized to its error o for a transverse momentum cutoff' of 1 GeV (open circles) and 2 GeV 
(closed triangles) of the D* . (b): Yield of the D* signal with p T (D*) >1 GeV (open circles) 
and >2 GeV (closed triangles) respectively relative to a momentum cutoff ofp T (D*) >2.5 GeV. 



3.2 Charm- Anticharm Correlations 

For many applications and QCD studies the two-particle inclusive distribution is of interest. 
We investigated the probability to fully reconstruct two charmed particles in HI. 

The acceptance that both charm quarks fall inside the acceptance of the tracking system 
(\r)\ < 1.5, p T > 1 GeV) is 3%. Using all the above decay channels and a vertex detector we 
can reach a double tagging efficiency of order 2.5xl0 -5 , which translates into 7000 charm pairs 
for 300 pb~ l . Without a vertex detector this number will be an order of magnitude less. 

If we do not require the full reconstruction of both charmed mesons, but ask only for one 
reconstructed D meson plus a lepton of momentum p T > 1.5 GeV which does not come from 
the primary vertex, we reach similar efficiencies as above, therefore doubling the number of 
double charm tags. Relaxing the quality of the double tagged events even more one would ask 
for one fully reconstructed charmed meson and a well separated jet. Efficiencies depend very 
much on chosen cuts and no numbers are given here. 



4 Summary and Comparison with other Machines 

At HERA, a rich physics program can be achieved with an integrated luminosity C = 300 pb^ 1 . 
The inclusive properties of open charm photo- and electroproduction can be studied at an ex- 
cellent level of accuracy. The next-to-leading order QCD predictions for the direct color-singlet 
photoproduction of J / ip can be thoroughly tested. Correlations between charm and anticharm, 
open bottom production, resolved and direct color-octet quarkonium photoproduction could be 
studied with £ = 300 pb~ l as well, although larger luminosity would allow a more significant 



no vertex 




acceptance 


reconstr. 


c -> £>* -> 




total 


detector 




p T > 2.5 GeV 


efficiency 


D U 7T 


decay mode 


efficiency 


D° -> K-7T+ 




0.03 


0.5 


0.18 


0.038 


l 10 4 












0.075 














0.032 




£)0 _^ 


-> 47T 








0.019 




sum (D°) 


0.03 


0.5 


0.18 


0.164 


4 10" 4 


with vertex 




acceptance 


reconstr. 


c^D°/D+ 


D°/D+ 


total 


detector 




p T > 1.0 GeV 
d/cr > 2.3 


efficiency 




decay mode 


efficiency 


D° -> ir-7T + 




0.05 


0.3 


0.55 


0.038 


3 10" 4 




~7T 








0.075 




D° -> 










0.032 






— > 47T 








0.019 




sum (-D ) 


0.05 


0.3 


0.55 


0.164 


l 10" 3 


D + -> - 


> 3tt 


0.05 


0.3 


0.21 


0.019 


6 10~ 5 




f 








0.091 






~7T~ — > 57T 








0.024 




sum (D*) 


0.05 


0.3 


0.21 


0.134 


4 10" 4 


sum (all D-mesons) 


2 10" 3 



Table 1: Summary of various reconstruction efficiencies for D° mesons in photoproduction 
with and without a vertex detector. Charged particles have been confined to \t)\ <1.5, and a 
vertex separation ofd/a >2.3 has been assumed for the lower example with a vertex detector. 
Acceptance and reconstruction efficiencies depend on details of cuts and only rough numbers 
using cuts similar to ref. are given. 



tagging method 


charm-tagging bottom-tagging 


D*-meson 

identified muon 

non- vertex track, djo > 3 

non- vertex track, djo > 2 


0.014 0.09 
0.073 0.20 
0.10 0.33 
0.37 0.17 



Table 2: Summary of charm and bottom tagging efficiency inside jets with E ] T et > 6 GeV, 
\r] jet \ < 2 , Q 2 <4 GeV 2 , and 135 GeV < W 1P < 270 GeV. Non-vertex track means a track 
with an impact parameter d larger than 2 resp. 3 times its error o (from ref. ffiffl). 



comparison between theory and experiments. Charm correlations can be exploited in order 
to measure the gluon density in the proton. In polarized ep collisions, charm data are useful 
to constrain the polarized gluon density in the proton. Very large luminosity is mandatory in 
order to study the production of \i J + 7 an d to investigate the charm fragmentation into 
J/ip. Intrinsic charm signals are not detectable without upgrading the detectors in the very 
forward region. A vertex detector is essential to reach a high tagging efficiency for open charm 
and should be taken into account in the design of the luminosity upgrade program. 



A figure of merit of the potential for charm physics is the total number of reconstructed open 
charmed particles. The biggest competition comes from fixed target experiments at Fermilab. 
A program to reach 10 6 reconstructed open charm is in place which corresponds to 300 pb^ 1 at 
HERA. Progress beyond 10 6 will depend on commitment to the physics and a potential to get 



10 at Fermilab is enticing |35| 
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